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ABSTRACT 


The  field  equation  for  meridionally  inextensional  axisymmetric 
deformations  of  shells  of  revolution  is  developed  and  applied  to  the 
postbuckling  analysis  of  simply-supported  imperfect  annular  plates  sub¬ 
jected  to  compressive  edge  loads.  The  semi-inextensional  model  is  then 
compared  with  the  usual  extensional  model  by  using  an  energy  approach 
with  two  and  three  deformation  parameters. 

The  method  is  used  to  analyse  the  postbuckling  behavior  of  an  imperfect 
annular  plate  with  an  eccentric  interior  edge  stiffener.  The  results  show 
that  although  either  imperfection  or  eccentricity  lead  to  stable  behavior, 
their  combination  may  introduce  unstable  imperfection-sensitive  behavior 
in  some  cases. 


TABLE  OF  CONTENTS 


Abs tract 

Table  of  Contents 

1.  Introduction 

2.  Plate  Model,  Geometry  and  Kinematics 

3.  Equilibrium  Equations 

4.  Constitutive  Laws 

5.  Approximations 

6.  Semi-Inextensional  Axisymmetric  Shell  and 
Plate  Analysis 

7.  Three-Parameter  Extensional  Model 

8.  Radially  Inextensional  Model  with  an  Eccentrically 

Reinforced  Edge 

9.  Buckling,  Post  Buckling  and  Sensitivity  of  Eccentrically 
Stiffened  Imperfect  Annular  Plates 

10.  Summary  of  Results 

Appendix:  Equations  for  a  Three-Parameter  Semi -Inextensional 
Model 

References 
List  of  Figures 


1.  INTRODUCTION 


The  postbuckling  deformations  of  many  plates  and  shells  exhibit  rather 
large  rotations  and  curvature  changes  while  the  midsurface  extensional 
strains  have  a  rather  limited  influence  on  the  deformation  patterns.  This 
type  of  behavior  is  related  to  the  so-called  "small-strain-finite-rotation" 
theory  of  plates  and  shells  which  assumes  that  the  equations  can  be  linearized 
in  the  extensional  strains  but  not  in  the  rotations  (and  displacements). 

In  many  practical  applications  it  has  been  observed  that  the  incremental 
extensional  strains  during  the  post-buckling  deformation  process  are 
sufficiently  small  and  over  most  of  the  shell  (or  plate)  so  that  the  de¬ 
formation  could  presumably  be  approximated  by  an  inextensional  model 
which  eliminates  the  incremental  extensions  completely.  However,  attempts 
to  do  so  have  been  found  to  be  not  entirely  satisfactory  because  of  two 
main  reasons: 

(a)  The  need  for  the  existence  of  localised  regions  of  high  strain 
gradients  in  order  to  account  for  the  boundary  conditions  and  for  the 
internal  compatibility  of  the  shell  deformation  pattern.  A  good 
example  is  the  diamond  pattern  of  the  buckled  very  thin  cylindrical 
shell  which  requires  significant  strain  gradients  along  its  creases  and 
corners  where  the  Gaussian  Curvature  of  the  surface  undergoes  large  changes. 

(b)  In  some  cases  the  inextensional  deformation  has  a  definite 
directional  pattern.  That  is,  the  incremental  deformations  are  virtually 
inextensional  in  one  direction  at  each  point,  while  having  a  substantial 
extensional  component  in  the  direction  orthogonal  to  it.  Quite  a  few 
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postbuckling  patterns  possess  this  property.  To  cite  a  few  examples,  the 
strong  postbuckling  of  thin  plates  in  shear,  of  thin  shells  of  revolution 
in  twist,  of  cylindrical  shells  and  of  annular  plates  in  axisymmetric 
compression,  are  all  of  this  type. 

For  the  latter  class  of  structures,  semi-inextensional  models  can  be 
postulated.  By  doing  so,  the  postbuckling  analysis  can  be  simplified  both 
by  decreasing  the  number  of  variables  and  by  a  priori  taking  into  account 
the  observed  deformation  patterns.  In  some  other  cases,  a  combination  of 
semi-inextensionality  plus  localised  areas  of  high  strain  gradients  may 
be  useful. 

Some  of  the  features  of  semi-inextensional  plate  models  will  be 
exemplified  in  this  study  by  considering  the  axisymmetric  postbuckling 
behavior  of  an  imperfect  annular  plate.  In  addition,  the  study  will 
include  an  examination  of  the  interaction  of  shape  imperfections  and 
boundary  eccentricities.  This  might  be  of  some  use  to  the  understanding 
of  this  complex  phenomenon  in  more  complicated  structures. 

The  problem  of  the  axisymmetric  postbuckling  behavior  of  annular  plates 
has  been  considered  in  the  past  by  several  investigators  [7,9,10,11.12,15] 
with  the  clamped  external  boundary  being  the  more  popular  case.  Treatment 
consisted  of  solving  the  axisymmetric  Von-Karman  plate  equations  by  using 
eith  series  expansion  techniques  or  especially  devised  numerical  methods. 

In  recentyears,  general  and  special  purpose  computer  codes  were  utilised 
too.  These  methods  have  their  advantages  in  accuracy  but  suffer  from 
the  drawback  of  being  difficult  to  use  for  overall  design  data  or  for  the 


evaluation  of  trends. 
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The  approximate  analytical  approach  prsented  here  for  the  evaluation 
of  simply  suported  annular  plates  under  compressive  edge  loads  is  based 
on  experimental  observations  both  for  its  basic  assumption  of  radial 
inextensionality  and  for  its  choice  of  deformation  model.  The  semi- 
inextensionality  assumption  is  used  for  deriving  the  differential  equa¬ 
tion  for  the  large  rotations  of  shells  and  annular  plates,  whereas  the 
deformation  model  is  combined  with  an  energy  approach  to  produce  a  simple 
approximate  model  of  post-buckling  behavior.  It  is  believed  that  this 
and  similar  approaches  are  suitable  for  engineering  estimates  and  as 
necessary  adjuncts  to  testing  programs. 

Observations  made  during  a  recent  testing  program  [8,p.l9;  10]  on 
simply  supported  annular  plates  have  demonstrated  that  as  high  post- 
buckling  loads  were  reached,  plates  with  free  internal  edges  or  with  weak 
rings  tended  to  deform  into  snapes  which  were  close  to  shallow  cones. 
Marked  deviations  appeared  only  if  the  stiffening  rings  at  the  interior 
edge  were  relatively  heavy.  A  similar  mode  of  behavior-  although  with  a 
different  loading  case  -  was  noted  in  the  study  of  shallow  conical  springs 
by  Almen  and  Laszlo  in  1936  and  later  by  Wempner  [2]  and  by  Schmidt  and 
Wempner  [3].  This  observation  will  be  subsequently  used  in  the  determina¬ 
tion  of  the  initial  deformation  model. 


s 
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2.  PLATE  MODEL,  GEOMETRY  AND  KINEMATICS 


Let  (r,z,6)  be  a  cylindrical  coordinate  system  in  space,  to  which  an 
annular  plate  or  shell  is  referred  (see  Fig.  1).  The  interior  boundary  of 
the  plate  (r=b)  may  be  stiffened  by  an  eccentric  ring.  Compressive  radial 
forces  P  (per  unit  lengthjare  applied  to  the  exterior  boundary  (r»a) . 

The  initial  configuration  of  the  imperfect  plate  is  that  of  a  very 
shallow  shell  of  revolution,  with  <J>Q(r)  being  the  angle  between  the 
normal  to  the  surface  and  the  positive  direction  of  the  z  axis.  (Note: 
the  magnitude  of  $  in  Fig.  1  is  exaggerated.)  Length  measure  along 
the  plate  midsuri.ace  in  the  r-z  plane  is  denoted  by  s  and  its  local 
radius  of  curvature  by  r^ .  Obviously,  the  following  two  relationships 
hold: 


1_ 
r . 


d<p 
_ c 

ds 


and  dr 


cos<|!  ds 
o 


(1) 


The  plate  deforms  from  its  initial  configuration  into  another  shell 
of  revolution  with  length  measure  s' .  The  radial  displacement  of  a 
material  point  on  the  plate  midsurface  is  denoted  by  u  (Fig.  2).  The 
angle  between  the  normal  fr  the  deformed  plate  and  the  z  axis  is  denoted 
by  $  ,  and  the  radius  of  curvature  of  the  deformed  plate  midsurface  is 
r ' .  The  following  relationships  hold: 


1  „  <&- 
7^  ds ' 


dr'  =  cos4>ds' 


Note  that  primes  relate  to  the  deformed  configuration  (except  for  d> 
which  is  measured  in  the  deformed  conf  igurationj , 
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An  expression  for  the  meridional  strain  can  be  derived  directly 
from  Fig.  2: 


ds'cos®  =  ds  coa<p  +  du 
o 


or, 


in  view  of 


e 

r 


(1): 


-  (1  +  t 


CO  sc 

o 

COSd 


-  1 


(2) 


See  also  equation  (15)  in  Ref.  [lj. 

The  radial  displacement  can  be  calculated  from  (2)  as  follows: 
r 

u  ®  f  [  (1  +  t  )cos*>  -  cos f  ]ds  +  u  (3) 

r  o  o 

) 

The  circumferential  strain  e  is  obtained  directly  from  its  definition 

c 

%  =  a/r  («) 

These  also  lead  to  the  strain  compatibility  equation: 

(1  +  ■  )cos(f  -  cosC  =  4—  (re.)  (5) 

r  o  ds  C 

Bending  strains  should  be  derived  in  principle  from  the  usual  Kirchoff- 
Love  assumptions  of  transversely  rigid  shells.  There  is  no  final  agreement 
on  the  "proper"  form  to  be  used.  Strict  application  of  the  K.L.  hypothesis 
leads  to  a  formula  of  the  type: 

K*-(l  +e)(|r  “  f“)  (6a) 

Here  e  is  an  extensional  strain  and  R  is  a  radius  of  curvature. 


On  the  other  hand,  an  alternate  definition 


-6- 


K  .  (1  «>£r  -  i 


(6b) 


is  widely  used.  It  has  also  been  shown  to  have  a  wide  range  of  validity 
in  small  and  large  strain  analysis  [14,  p.  94],  Their  difference  is  the 
quantity  e/R.  Fortunately,  in  this  study  the  strains  are  small  and, 
moreover,  the  imperfect  undeformed  plate  is  almost  flat  so  that  R  is 
very  large.  Hence,  the  difference  is  quite  negligible,  and  the  second 
definition  will  be  adopted  for  convenience.  This  leads  to  the  following 
expressions  for  the  bending  strains: 


K 

r 


ds '  ,d$ 


d<P 


ds  ds 


(££.  ) - 2.  =  ) 


ds 


Ke 


r+u  ^siny ^ 
r  v  r+u 


sin$ 
_ o 

r 


— (sin<j>-sin<j)  ) 
r  o 


(7) 

(8) 


This  is  also  the  form  adopted  (as  an  approximation!)  by  Reissner  in  his 
nonlinear  theory  of  shells  of  revolution  [1,  eqns.  17,18], 

From  equations  (3),  (4),  (7),  (8)  it  follows  that  all  the  extensional 
and  bending  strains  have  been  expressed  in  terms  Ef,  <(>  and  the  constant 
of  integration  uq. 


3.  EQUILIBRIUM  EQUATIONS 

Equilibrium  is  considered  in  the  deformed  (prime)  configuration.  The 
force  resultant  vector  acting  on  the  cross  section  r'  =  c  can  be  cecomposed 
into  either  tangential  (NM  and  transverse  (Q')  or  radial  (H' )  and  axial  (V') 
resultants  respectively.  These  are  related  by  the  equations  (see  Fig.  2): 
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V*  -  N^sin$  +  Q'cos<ti 


H'  =  N^cos4>  -  Qrsin$ 


Axial  equilibrium  requires  that  V'  = 
terms  of  Q*  as  follows: 


0,  so  that  can  be  expressed  in 


N’  =  -Q’cot$ 
r 


The  second  equilibrium  equation  is  in  the  radial  direction  and  yields 
the  expression: 


zr(HV>  -*e 

Now,  H'  can  be  eliminated  from  the  equations  by  using  the  expression 
for  and  substituting: 

H'  =  -  0  ’cot<?cos<f>  -  Q’sinif  =  - 

S  i.Il(P 

go  that  the  equilibrium  equation  becomes 

srO +  %  - °  (i2> 

In  many  applications  stress  resultants  defined  per  units  of  length  of 

the  undeformed  shell  cross  sections  are  of  advantage,  especially  in  conjuntion 

with  constitutive  laws.  These  are  defined  by: 

Q  =  —  Q'  ;  N  -  —  N'  ;  M  -  —  M’  ;  N  «  N'  ;  M.  =  ~~  K 
r  r  r  r  r  r  r  8  ds  9  0  ds  6 


Then  the  force  equilibrium  equation  reduces  to 


6.  (Qf-)  +  N  =  0 

ds  sin<j>  9 


(13) 


The  equation  of  moment  equilibrium  in  the  prime  configuration  is  that  of 
a  shell  of  revolution: 


4r  (r'M1  )  -  M’r'cos-J!  -  r'r'Q' 
d<p  r  d  <j>  (j> 

or: 

-  M’cos»  -  r'Q' 

Passing  again  to  the  stress  :esultants  per  unit  undeformed  length, 
the  following  results: 


--U-M  )  - 
as  r 


M  cost  - 


ds' 

ds 


rO  = 


(1+t  ) rQ 
r 


(14) 


Hence,  all  the  equations  ot  equilibrium  are  reducible  to  the  single 
equation: 


.  ! I M  >  -  fosiil  ; 

C_  ills _ t _ 9  ! 

ds  !  (1  +  ?.  )smj  !  0 

r  -J 

Note  that  all  the  variaoj.es  ate  measured  now  with  respect  to 
of  the  undelonned  configuration. 


(15) 


the  geometry 


4 .  CONSTITUTIVE  LAWS 

It  will  be  assumed  that  the  stress  resultants  per  unit  undeformed  length 
are  related  to  the  strain  measures  by  the  usual  Hooke’s  law.  This,  plus 
the  other  assumptions  of  the  Kir chho f f -Love  hypothesis  lead  to  the 
commonly  used  constitutive  laws  of  plate  and  shell  theory: 


-9- 


M  =  -D(K  +  vKfl ) 
r  r  6 

Mq  =  -D(K0  +  vKr) 

Nr  =  C(£r  +  V£e) 

Ne  -  C(eQ  +  vef) 


3W 

3K6 

3W 

3c 


(16) 


with 

C  =  — —  :  D  =  -ai-  -  £  C 

1-'/  12 (1-v2 )  12 

The  strain  energy  function,  per  unit  area  of  the  undeformed  reference 
surface  is  also  given  in  its  common  separated  form 


W 


+•  c2+ 


t  +  Ke  +  2vKrKe> 


(17) 


Here  the  shear  terms  have  been  deleted  because  of  the  axisymmetry.  Re¬ 
garding  the  magnitudes  of  the  strains,  equation  (17)  should  be  valid  for 
the  same  range  of  (small)  strains  for  which  Hooke's  Law  is  valid.  For 
larger  strains,  more  refined  forms  for  the  strain  energy  (or  else,  other 
constitutive  laws)  ought  to  be  used. 


5.  APPROXIMATIONS  (GENERAL  CONSIDERATIONS) 

There  are  several  approximations  and  simplifications  to  the  nonlinear 
expressions  obtained  above.  These  are  related,  as  is  usual  in  plate  and 
shell  theory,  to  four  distinct  properties  of  any  given  shell  problem: 


* 
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Ca)  Constitutive  laws  and  the  magnitude  of  the  strains.  Here,  the  small 
strain  approximation,  which  is  usually  made,  is  valid  for  most  engineering 
materials  in  the  elastic  (and  low  plastic)  range. 

(b)  Particular  geometry.  Here,  it  may  be  assumed  that  the  initial  imperfec¬ 
tion  <P  is  sufficientlv  small  so  that  the  approximation  dr  %  ds  can  be 

o 

(carefully!)  introduced  where  appropriate. 

(c)  Restrictions  on  the  magnitude  of  the  rotation.  Here,  the  quadratic 
approximation  can  be  introduced  for  Cos>)>  (to  within  0.4%  for  <)>  <_  30°)  and 

the  linear  (to  within  4.5%  for  4>  ^  30°)  or  cubic  (to  within  .07%  for  <p  <_  30°). 
approximations  introduced  for  sin<p. 

(d)  Assumptions  based  on  the  estimated  mode  of  behavior.  Here,  the  semi- 
inextensional  approximation  may  be  introduced  by  taking  e  -*■  0.  This 
assumption  must  be  justified  by  apastaaiAti  comparisons  with  more  complete 
solutions  oc  with  experimental  observations.  It  should  be  noted  here,  though, 
that  this  is  not  the  only  possible  approximation.  For  example,  the  methods 

of  references  four  and  thriteen  which  omit  or  replace  terms  in  the  elastic 
strain  energy  function  have  been  used  with  a  varying  amount  of  success  for 
the  study  of  thelarge  deflection  of  plates. 

The  present  study  makes  selective  use  of  the  first  three  approximations, 
with  the  intention  of  developing  and  examining  the  fourth,  e.g.  the  semi- 
inextensional  assumption. 

Some  examples  of  the  uses  of  approximations  (a),  (b)  and  (c)  are  as 


follows: 
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Using  assumption  (b)  and  (c),  the  expression  for  u  becomes 

u  =  f  [  (l+e^)cos<{>  -  cos$0]dr  =  J[(  1-  j4>2)Er  +  j(4>^-4>2)  ]dr  +  0(<)i4) 


The  expression  for  in  terms  of  u  and  $  is: 

e  •  (1  +  4^  )cos4)  sec^  -  1  =  4^  +  4^‘f’2-^2)  ••• 

r  dr  o  dr  z  o 


whereas  the  compatibility  equation  reduces  to: 


4“(rO  =(1  -  +  j($2  -b2)  +  0(<fr4) 

dr  0  2  r  2  o 


Using  (a)  and  (b),  the  expression  for  the  meridional  bending  strain  is: 


K  =  ~(>P-0  ) 
r  dr  o 


Using  (a)  and  (b),  the  curvature  compatibility  equation  reduces  to: 


4-<iM  ~  cosoK 
dr  0  r 


Introduction  of  the  quadratic  approximation  (c)  for  cos<j>  yields  further: 


3prV  ‘  a  '  r»X  +  °‘*5) 


This  equation  is  of  the  same  order  of  accuracy  as  the  cubic  approximation  for 
K  ,  together  with  (b): 


cKg  -  (♦  -  <Po)  -  |  (*3  -  4>^)  +  0(e) 


For  most  purposes  the  linear  approximation  is  sufficient: 


rK0  =  4>  -  <t>o  +  0  (4>  ) 


so  that: 


f?(rKe>  "  Kr 


6.  SEMI  -  f.  NEXTEN  S IONAL  AXISYMMETRIC  SHELL  AND  PLATE  ANALYSIS 


In  this  model,  the  radial  strains  are  omitted  in  the  shell  analysis 
in  comparison  with  the  circumferential  strains  and  the  rotations  (whose 
magnitudes  are  not  restricted). 

Dif ferentiation  of  the  equilibrium  equation  (15)  while  making  use  of 
(2),  (4)  and  (16),  and  pucting  e  «  0  throughout,  yields  the  expression: 


1  f  d  \h  UMr 

is  "^dsT  [  s 


-r~  (rM  )  -  cosJMfl') 

— x - -Hr  +  C(cos?  -  cos'?  )  -  0 


Making  further  use  of  (16)  and  also  of  (7)  and  (8)  leads  to  a  fourth  order 
differential  equation  in 

d  i  f  d  f  i'd,  d?  d$ox  cos$  ,  ,  A  A  x  , 

ds  ,  I.  ds  I  sin  L.  ds  ds  ds  r  o 

,  •j’l  ||  ^2 

+V  - —  (tosf  -  cos?>  )  L  -  (cos<?  -  cos<?  )  =  0 
ds  of,,,?  o 


Noting  that  for  thin  plates  and  shells  —  >>1,  the  omission  of  the  term 
d<t>  h2 

with  in  comparison  with  the  last  term  is  justified,  yielding  the 

simplified  form: 


^  (sln *  -  sinV}jr  j(cos^o8^0>  =  o 

a  L  L  )  , )  h 


(18) 
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The  last  equation  describes  the  axisymmetric  serai-inextensional 
deformations  of  shells  of  revolution  in  general.  The  input  into  the  equa¬ 
tion  is  the  initial  shape  [r(s)  and  $  Cs ) ]  and  appropriate  boundary 

o 

conditions.  The  result  would  be  $(s)  which  expresses  the  angular 
coordinate  of  the  deformed  shell  interms  of  the  Lagrangian  variable  s. 

The  other  terms  which  complete  the  description  of  the  deformed  configuration 
are 


u  =  f (cos?  -  cos<{>  )ds 
'  o 

v  =  |(sin(t>  -  sin4>  )ds 
o 

where  u  and  v  ate  the  radial  and  axial  displacements  respectively. 

The  bending  strains  (K  ;K  )  and  extensional  strain  eu  can  then  be  calculated 

v.  y 

from  (7),  (8)  and  (4),  thus  facilitating  the  calculation  ofthe  moments 

(M  ;  M  ),  force  resultant  N  ,  and  shearing  force  Q  .  The  meridional  stress 
r  “  r  r 

resultant  N  can  be  calculated  from: 
r 

N  =  -  Q  Co  tip 
r  r 

(The  constitutive  laws  cannot  be  applied  for  the  calculation  of  because 

of  the  meridional  inextensionality  assumption).  Note  that  for  $  *  0, 

cannot  be  calculated  in  a  direct  fashion  because  of  the  apparent  singularity, 

but  lim  N  does  exist,  as  will  be  seen  later 
r 

4>-*o 

The  imperfect  annular  plate  is  included  in  the  broader  class  of  "shallow 
initial  configurations"  defined  by  approximation  (b).  Replacing  ds  with  dr, 


the  result  is: 


-2.  (Tr  - 2.  [r  _a  *  )] 

dr  H  dr  ]  sin<t>  dr  1  dr  KV  VJ 


-  (aint  -  sin4>  'L  }}  - 

r  o  • 


12 


2  (cosij>  -  cossj>o)  =  0 


(19) 


Here,  also,  the  following  can  be  used  (if  desired) : 


\ 


sini}>  y  d> 
o  '  o 


$ps*0  % 1  -  j  i  % 1 


The  latter  approximations  do  not  lead  to  any  obvious  advantage  unless  the 
deformed  configuration  is  shallow  too,  in  which  case  an  expansion  of  the 
trigonometric  terms  up  to  quadratic  (or  cubic)  terms  in  <f>  would  make  the 
equation  more  "manageable"  from  a  mathematical  point  of  view.  In  such  a 
case,  the  use  of  a  new  variable 

3  =  $ 

o 


and  expansions  in  8  may  be  of  advantage. 

Equations  (18)  and  (19)  can  be  cast  into  algebraic  forms  by  using  the 

new  variable  y  =  sin®. 

1 

.  n  2.2 

cosif)  =  (1  -  y  ) 


i!  = 

A  a 


a 


i 

2,2  d£ 

y  )  if 


Then: 
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and  so  forth.  This  form  is  particularly  suited  to  uses  of  the  equations 
in  conjunction  with  standard  numerical  programs  for  solving  ordinary 
differential  equations. 

Methodology: 

While  further  analysis  of  the  semi-inextensional  differential 
equation  may  be  rewarding  in  itself, the  purpose  of  comparing  this  model 
with  the  more  complete  (and  more  complex)  extensional  model,  can  be 
achieved  by  far  simpler  means  which  should  offer  their  own  rewards. 

At  the  first  stage  of  the  investigation,  simle  deformation  patterns 
will  be  used  in  order  to  compare  the  regular  plate  model  (which  includes 
the  effects  of  e  )  with  the  semi-inextensional  model.  The  theorem  of 

L 

the  minimum  of  the  total  potential  will  be  utilized  to  obtain  the 
response  of  the  plate  to  compressive  edge  forces  within  the  assumed  de¬ 
formation  modes.  At  this  stage  the  inner  edge  of  the  plate  will  be 
free.  Hre  initial  imperfection  of  the  plate  will  be  in  the  form  of  a 
very  shallow  cone  with  a  constant  .  The  fact  that  this  is  actually 
the  governing  imperfection  of  the  simply  supported  annular  plate  has  been 
borne  out  by  experiments  [8,  Figs.  9b,  17a,  17b], 

The  comparison  between  the  two  models  should  provide  an  indication 
of  the  adequacy  of  the  semi-inextensional  model,  since  a  three  parameter 
extensional  model  will  be  compared  with  a  two  parameter  semi-inextensional 
model.  The  comparison  will  concentrate  on  the  region  of  the  buckling  point 
itself.  The  reason  is  that  the  semi-inextensional  model  should  have  its 
greatest  difficulties  and  inaccuracies  as  4>,4>  -*0. 
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A  three  parameter  more  complex  semi-inextensional  model  will  alao  be 
presented  (see  Appendix).  It  will  not  be  investigated  in  detail  since  it 
has  been  argued  that  only  the  simplest  nontrivial  models  should  be 
investigated  by  analytical  tools  and  more  complex  models  are,  in  many 
cases,  inferior  to  direct  numerical  procedures.  This  point  is  left  open. 

At  the  next  stage,  the  effects  of  an  interior  edge  beam  will  be  incorporated 
into  the  semi-inextensional  analysis.  This  should  facilitate  the  study  of 
the  comparative  roles  and  interactions  of  the  edge  beam  eccentricity  e 
and  the  initial  imperfection  4>  .  The  possibility  that  the  interaction 
can  cause  essential  changes  in  the  response  will  be  looked  into.  The 
fact  that  an  eccentrically  stiffened  imperfect  plate  can  be  imperfection 
sensitive  has  been  verified  experimentally. 


7 .  Three-Parameter  Extensional  Model 


Consider  a  displacement  field  of  the  form: 

u  =  u  t  ru, 
o  1 


(20) 


where  is  the  initial  Imperfection  and  uq;  u^;  are  constant  de¬ 
formation  parameters  to  be  determined  by  minimizing  the  total  potential. 

The  deformation  pattern  assumes  that  the  dominant  post-buckling  mode  is  that 
of  a  shallow  cone.  This  has  been  well-correlated  with  experiments  on  annular 
plates  [8,  p.  19  and  figs.  9b,  17a,  17b].  Using  the  quadratic  approximations 
for  the  various  strain  measures  and  applying  equations  (2a),  (4),  (7),  (8b), 
the  following  expressions  result: 
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er  “  U1  +  b*l*P 


Also: 


-  u  +  u, 
r  o  1 


K  =  0 

r 


Kr  =  — (?.  -<P  ) 
0  r  1  o 


The  total  elastic  potential  of  the  system  is: 
a  a 

V  -  nof  (K2+K.2+2vK2K.2)rdr  +  ircj(e2+e2+2ve  ea)rdr  +  2iraPu  (21) 

S  r  6  r  0  •  r  ©  r0  a 

b  b 

where  u  =  u(r=a)  and  P  is  positive  in  compression. 

cL 

introducing  the  values  of  tg;  K^;  K^,  the  result  for  V  is: 

V  =  rrD(<j>  -®  )2in(a/b)+  ncj^(a2-b2 )  [u  +  ■5L(<j»2-4>2)  ]2  +  [u2ln(a/b)  + 

lo  /  4  1  2  1  o  o 

+  2u  u  (a-b)  +  ~u2(a2-b2)J  +  2v{u.  +  ]  [u  (a-b)  +  \  u..(a2-b2)}} 

Oi  21  12100  41 

+  2iraP  (u  +au, )  (22) 

o  1 


The  total  potential  can  be  now  minimized  with  respect  to  the 

parameters  u  u, ,  <p  ,  leading  to  three  simultaneous  equations  in  the 

o,  1  1 

parameters.  These  can  be  brought, after  some  manipulations  and  an  elimina¬ 
tion  process,  to  the  following  cubic  equation  for  : 


-  g  "  g2**i~*P*i 


(23) 


where  the  loading  parameter  \  and  modified  rotations  <t>^  are  given  by: 


Pa2 

D 


*i  ‘  h  *i 


i 
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Also: 


g 


1 


-  2  [  (1+y  )  £nv+ ( 1+v )  ( 1-ti )  ]  fcny 
1-U  2+v ( 1-y ) 2+(l+u ) i ny 


g2^> 


— 6(l-u)  (1-Vi2)  [1-P+  -j(l+y)£.ny  ] 
l-y2+v(l-y )2+(l+y)&nu 


y  =  b/a 


For  a  given  loading  X,  geometrical  ratio  y  and  imperfection  <t>o> 
the  equation  can  be  easily  solved  for  Once  is  known,  the  displace¬ 

ment  parameters  u^  and  can  be  easily  calculated.  These  are  given  by: 

uo  -  g3(^-r>  +  g4A  (24) 

U1  *  +  r6A  <25> 


where: 

g3  -  ij  (1-v)(1-ij2) 

-  ^ 
g4  =  "  124 

85  *  f^fvU-u)  +  \  (1+uHny] 

1  r .  .  In y  , 

86  =  124  [1  (1-y) (1+v)  ] 

A  -  -  (1+yHny  -  (1+v)  (1-u) 


u 


i 


h2 


u 


i 
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The  special  case  <(>  =  0  is  that  of  the  perfect  plate.  The  equation 

for  <f^  breaks  down  into  the  straight  line  =  0  which  describes  the 
prebuckling  behavior  and  the  parabolic  "equilibrium  path": 

X  =  h  ~  *2*1 

which  describes  the  postbuckling  behavior  of  the  plate.  The  intersection 
is  at 


where  is  the  critical  (buckling)  loading  parameter.  The  equilibrium 
paths  of  the  imperfect  plate  are  cubic  parabolas.  These  parabolas  approach 
the  two  curve  system  of  the  perfect  piate  as  <Pq  *  0.  Also,  for  a  given 
4^,  the  two  curves  approach  each  other  for  4>^>>4>o.  This  signifies  the  fact 
that  for  very  large  displacements  the  initial  imperfection  is  not  too 
important . 

Considering  now  the  radial  displacements,  it  is  convenient  to  refer 
to  the  displacement  at  a  given  point.  The  inner  edge  displacement  u^  -u  ( r=b ) 
is  chosen  because  of  the  direct  connection  between  this  displacement  and 
the  circumferential  stress  level  at  the  inner  edge.  An  expression  for 
is: 

^  =  (g3  +  wg5)  +  (g^+uggH  (26) 

Another  expression  can  be  obtained  by  eliminating  through 

the  use  of  the  equation  of  the  cubic  parabola.  This  yields: 
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The  equilibrium  path  of  the  perfect  plate  $  =  0  consists  of  two  straight 

lines:  the  prebuckling  line 


“b  =  (g4+vl86)A 


and  the  postbuckling  line 

=  ij  ~  -  84  “U86)X 

The  two  lines  intersect  at  the  buckling  point,  where 


^  V  =  (g4+^6)8l 

Note  that  the  radial  motion  reverses  its  direction  at  the  buckling 
point.  It  is  inward  up  to  the  point  and  outward  thereafter. 

Again,  the  equilibrium  path  of  the  imperfect  plate  provides  a  smooth 
transition  between  the  pre-  and  post-buckling  regions,  approaching  the 
straight  line  of  the  perfect  plate  for  . 

The  fact  that  the  high  postbuckling  curve  of  the  plate  approaches  a 
straight  Line  was  observed  experimentally  [8],  [10],  where  the  use  of  the 
two  asymptotic  straight  lines  for  the  experimental  determination  of  the 
critical  load  led  to  satisfactory  results  [10],  This  agrees  quite  well 
with  the  results  of  this  section  and  indicates  that  the  three  parameter 
model  is  sufficient  to  indicate  the  main  features  of  the  deformation 
process.  As  will  be  seen  later,  the  model,  gives  also  surprisingly  good 
quantitative  predictions  of  the  buckling  point  over  a  wide  range  of  the 
geometrical  parameter. 
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8.  Radially  Inextensional  Model  with  an 
Eccentrically  Reinforced  Edge 


The  condition  of  radial  inextensionality  imposes  a  relationship 
between  the  rotations  and  radial  displacement.  For  a  quadratic  approx¬ 
imation  of  the  cosine  function  the  relationship  is: 


,  r 

u  =  ^  -  jf  (<pz-<t>z)dr 


(28) 


Hence: 

ce  '  X  -  tX-«o>dr1  <29> 

b 

The  expressions  for  K.  and  K.0  in  terms  of  the  rotations  remain  as 
before  (7a,  8b). 

It  is  assumed  that  an  eccentric  ring  with  area  A  moment  of  inertia 
I  and  eccentricity  e  (with  respect  to  the  plate  center  line)  is  attached 
to  the  inner  boundary  of  the  plate  at  r=b.  Positive  eccentricity  is 
defined  such  that  positive  rotations  at  r=b  would  cause  outward  motion 
of  the  extreme  fibers  of  the  ring  (Fig.  1).  Under  these  conditions, 
the  circumferential  strain  In  the  ring  is  given  by: 


'0R 


7  *  h\  + 


(30) 


where  z  is  the  axial  coordinate  in  the  ring  (positive  in  the  e  direction). 
The  radial  variations  within  the  ring  are  omitted  in  the  analysis,  for 
simplicity. 

The  inclusion  of  the  ring  in  the  analysis,  while  not  important  for  the 
comparison  with  the  extensional  model,  offers  the  capability  for  investigating 
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the  interaction  of  boundary  eccentr icities  with  geometrical  imperfections 
and  their  combined  effect  on  the  post  buckling  behavior. 

The  total  potential  of  the  system  of  a  radially  inextensional  annular 
ring  with  an  eccentric  stiffening  ring  along  the  interior  boundary  of  the 
plate  becomes: 

a  a 

V  =  ttC  J  rt2  dr  +  nD  1  (K2-HC2+2vK  K.)rdr  +  nEb/fc2  dA  +2Prtau 
b  0  b  6  r  r  6  A  0R  3 

=  r.C  J  jr  [ub~  j  '  (c‘'-r  )dr]2cLrtnD j  f  )  1 2  +  —  2  (^_<^0 > 2  + 


+  -  i-l  (0-^)2]  irdr  +  “  tKc-oo)2  +  IV^oV1^1 


+  27tal’lub  -  -y  / (<p‘  -$“)dr )  (31) 

b 

Here  A  is  the  cross  sectional  area  of  the  ring,  and  I  is  the  moment 
of  inertia  \  i  res.  et  t  to  a  radial  axis  through  its  centroid, 
sequel,  <f>  will  be  taken  to  be  that  of  a  shallow  cone.  This  form  of 
imperfection,  as  noted  before,  is  common  and  has  been  oh  served  experi¬ 
mental^. 

(a)  Two  Parameter  Model  with  a  Free  Edge 

Tnis  is  the  simplest  model  possible.  The  deformed  shape  is 
also  taken  to  be  that  of  a  shallow  cone.  It  is  similar  to  the 
three  parameter  model  discussed  above,  but  needs  only  two 
parameters  for  its  description.  The  parameters  are: 


Q  -  =  constant  and  u^. 
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Also,  I  »  A  =  0  are  taken  in  the  analysis  for  the  free  edge.  The  form 
of  the  strain  energy  simplifies  to: 

a 
f 
; 
b 


V  =  nC  J  — [u.  -  tr(r-b )  ($2-fi2 )  ]  2dr  4-  ttD  J  — ($  -$  )2dr  + 
^rbz  no  trio 

+  2.faP[ub  =  ^(b-a)(?i2-$2)] 


(33) 


or 


ita2C 


u  u,  . 

-  -in^C—-)2  -  (1-p+M inu) (— -  )  ($2 -4>?-)  +  -5-(3u2-^u+l-2iJ2inp )  C<t> ^  — <fs 2 ) 2 

a  3X00  X  O 

-  —  ini.u  -;  )2  +  -Al  A[^_  -  iu-w)(*?-*2)J  (34) 

12a2  1  °  6a2  a  2  1  0 


The  total  potential  can  now  be  minimized  with  respect  to  (u^/a)  and 
This  results  in  the  following  equations  for  ^  and  u^: 


'*1  =  W* o]  ~  f2al'^)<PL 

“b  -  £3(n4o}  +  v 


(35) 

(36) 


where: 


£,  = 


(in,.)  ‘ 


1  s- — 1— inu 

6(l-u ) :+  3(l-u? ) inn 


f2 


£ 


vi-l-tnu 
t-l-u int 


3  2tnn 

=  1 

4  12  ini- 
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An  alternate  expression  for  is  obtained  by  eliminating  : 


f,f.  f, 

-P  <-H-  >  -  P  -  £4)X 
2  V1  2 


The  similarity  in  form  between  the  3  parameter  extensions!  model 
and  2  parameter  semi-inexten9ional  model  is  self  evident.  Also  obvious  is 
the  fact  that  the  semi-inextens tonal  model  requires  a  smaller  number  of 
parameters  for  the  same  (rotational)  mode  shape  than  the  extensional  model. 
The  rest  of  the  analysis  follows  along  the  same  lines  of  the  3  parameter 
model,  with  the  "f"  functions  replacing  the  "g"  functions  and  with  g,_  and 
deleted.  The  buckling  load  is  given  by 


The  ptebuckling  line  roc  the  radial  displacements  in  the  perfect 
plate  is 


and  that  of  the  postbuckling  line  is: 


hh  h 

~2  "  l,Z  '  '4 


In  Table  1,  values  of  the  buckling  parameter  are  given  using: 

(a)  theoretical  (exact)  calculations  (b )  three-parameter  extensional  model 
and  (c)  two-parameter  semi-inextensional  model.  Also  given  are  the  other 
coefficients  of  the  postbuckling  equilibrium  paths. 


-25- 


Table  1:  Comparative  Values  of  A  and  Coefficients 


y 

0.2 

0.3 

0.4 

0.5 

0.6 

(a)  \heor 

3.55 

3.11 

2.77 

2.45 

2.30 

<b)Xext-*l 

3.68 

3.10 

2.75 

2.51 

2.35 

(c) 1  *f 

K  J  inext  1 

3,20 

2.88 

2.65 

2.49 

2.35 

(b)  g2 

m 

co 

i 

-.621 

-.436 

-.291 

-.180 

Cc)  f2 

2 

-.982 

_ 

-.688 

-.471 

-.309 

-.187 

i 

The  table  shows  that  the  extensional  model  gives  somewhat  better 

predictions  for  the  buckling  parameters  for  small  y.  This  difference 

disappears  for  increasing  y  and  both  models  yield  very  good  predictions. 

Even  for  y  =  0.2,  which  is  a  rather  extreme  case,  P.  is  within  10%  of 

mext 

the  exact  value. 

From  a  theoretical  point  of  view,  the  inextensional  model  is  not 
expected  to  yield  accurate  results  for  d>-*0  since  the  in-plane  extensions 
are  dominant  in  this  range  of  the  deformation  process  -  and  these  are  not 
accurately  described  by  the  model.  Accuracy  is  expected  to  improve  with 
increasing  <p.  Hence,  the  good  results  near  <j>=  0  are  quite  satisfying  and 
are  good  indications  of  the  overall  capability  of  the  model  even  at  extreme 
cases. 

(b)  Two  Parameter  Model  with  an  Eccentric  Stiffening  Ring 

The  relative  success  of  the  simpler  semi-inextensional  model  in  the 
case  of  the  free  inner  boundary,  makes  it  the  model  of  choice  for  the  more 
complex  case  of  the  eccentrically  etif fenedplate,  where  both  imperfections 
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and  eccentricity  interact,  it  should  be  stated,  though,  that  the  use  of 
a  two  parameter  model  restricts  the  analysis  to  relatively  weak  rings. 

The  reason  is  that  the  ring  provides  some  restraint  against  rotation  for 
the  inner  edge,  The  rotational  restraint  of  a  heavy  ring  is  sufficiently 
close  to  the  clamped  edge  condition,  so  that  the  "cone  approximation"  can¬ 
not  adequately  represent  the  bent  shape  of  the  plate,  and  a  minimum  of 
three  parameters  would  be  needed.  For  the  basic  equations  for  this  more 
complex  analysis,  see  Appendix. 

Introduction  of  the  assumed  deformation  mode  into  the  more  complete 
form  ot  the  total  potent  Lai  and  diiierentiaiion  with  respect  to  and  u^, 
lead  t.o  the  following  equations: 


(bl-40[d3+d6e  i-(3?,+f^ed2  i-  =  (e+^d^ 


(38) 


%  '  12CVV®  *  Ur* o)d2“Aj 


(39) 


where: 


A(l-,z) 

bh 


—  Ae 
e  -  -- 


d^  =  A  -  >nu 

d^  =  6(l-jj+v.enu) 

d^  *  (I  -  £nu)d^ 

d^  =  3 ( 3p 1 -4u+l-2u 2 £ n^ ) d^  -  ^ 
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dg  “  (l-lO(A-l)  -  Inn 

,  _  12£ny 

d6  -~T 

f.S. 

1  Db 

If  A  =  I  =0  are  introduced  here,  the  equations  for  the  free  edge 
result.  The  effects  of  the  ring  are  twofold:  Increases  in  the  area  A 
and  moment  or  inertia  1  of  the  ring  cause  a  corresponding  increase  in 
the  effective  stirrness  of  the  reinforced  plate  through  their  effects 
on  the  coefficients,  especially,  and  d^.  However,  these  do  not  lend 
to  qualitative  changes  in  the  behavior  of  the  plate.  On  the  other  hand, 
the  addition  of  the  eccentricity  leads  to  the  appearance  of  new  terms 
in  the  equations.  The  combination  of  eccentricity  and  initial  imper¬ 
fection  introduces  new  modes  of.  behavior  which  did  not  exist  inthe 
centric  plate.  This  will  be  studied  in  more  detail  in  the  sequel. 

8 ■  BUCKLING,  POST  BUCKLING  AND  SENSITIVITY  OF  ECCENTRICALLY 
STIFFENED,  IMPERFECT  ANNULAR  PLATES. 

The  investigation  is  carried  out  based  on  the  semi-inextensional  model 
with  two  parameters,  assuming  axisymmetric  behavior. 

As  is  well  known,  the  perfect  annular  plate  has  a  stable  initial  post- 
buckling  behavior  and  is  not  imperfection  sensitive.  The  addition  of  either 
an  eccentric  stiffening  ring  or  an  initial  imperfection  does  not  change 
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the  stability  of  the  behavior  of  the  plate.  The  structural  problem, 
though.  Is  changed  in  either  case  from  a  bifurcation  problem  plus  postbuck- 
ling  to  that  of  a  nonlinear  behavior  problem  with  the  external  compressive 
load  increasing  monotonously  with  the  deformation. 

However,  the  situation  is  changed  when  both  imperfection  and  eccen¬ 
tricity  occur  simultaneously.  It  can  be  shown  that  for  some  values  of  the 
imperfection  and  eccentricity  parameters,  the  structure  can  buckle  (that  is, 
it  has  a  bifurcation  point).  Furthermore,  this  point  is  unstable  and  is 
imperfection  sensitive  for  a  limited  value  of  the  imperfections  to  one  side 
of  it.  This  mode  of  behavior  should  be  anticipated  if  it  is  recalled  that 
the  eccentricity  reduces  the  order  of  Lhe  nonlinearity  from  the  cubic  to 
the  quadratic  [11]  and  that  structures  with  quadratic  nonlinearities  are 
usually  imperfec tion-sensicive  (the  different  behavior  of  frames  from  that 
of  beams  is  a  common  example). 

To  examine  this  mode  of  behavior  within  the  scope  of  the  simplified 
model,  it  is  convenient  to  rewrite  the  equations  by  introducing  the 
variable  and  "equivalent  Imperfection"  as  follows: 


fi  =  -  4> 

i  o 


1  - 
=  (f  +  m—  e 
o  d^ 


(40) 


The  equation  for  the  quilibiium  path  A  =  A (6)  becomes  then: 


X(8+e)  =  [d^6‘  -t-3(d ^e  +  d^c  )6  +  (dge2  +  4d^ec  +  2d^e2  +  d^) 


(41) 


with: 


d,  -  d,  -  -r-  «  3(1  +  2uinM  -u2)^ 
7  2  d5  d5 


d6  -h 


(d2  +  d7) 
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Note  that  for  0-  u  '  1,  all  the  coefficients  (d^...dg)  are  greater 
than  zero. 

When  the  "equivalent  Imperfection"  c  vanishes,  the  equation  reduces 
to  a  bifurcation  problem  with  the  two  branches 

3=0  (nreb  tickling) 

and 

a  =  (d£  +  3d^e3  +  dfie‘  +  d^)  (postbuckling)  (42) 

The  "buckling  loud"  ' 1  is  the  intersection  of  the  two  lines  and 

cr 

is  given  by: 


a’  =  ±-ld  ,  +  dQe->  (43) 

cr  d,.  ’  8 

The  formula  holds  also  fur  e  =  0  when  it  represents  the  "true"  buckling 
load  ■.  of  the  corresponding  centrically  stiffened  plate.  It  should 
be  stressed  again,  though,  that  its  validity  is  restricted  to  weakly 
stiffened  plates-  When  the  stiffener  is  relatively  heavy,  an  additional 


term  in  the  rperesentat  Lon  for  is  needed  and  the  two  parameter  model 


would  tend  to  exaggerate  the  buckling  load. 

When  e  i  0,  the  postbuckling  parabola  crosses  the  8=0  axis  at 
a  positive  slope.  This  implies  that  the  buckling  point  is  unstable  and 
the  load  decreases  foL  8  0.  The  minimum  of  the  parabola  is  reached  at 

9d7 

=  A«  -  - - —  g'' 

min  cr  4d,dc 
4  5 


This  value  of  \  .  is  higher  than  the  critical  loading  X  of  the 
min  °  cr 

corresponding  centric  plate,  indicating  that  for  small  eccentricities 
the  degree  of  sensitivity  of  the  postbuckling  behavior  is  not  very  severe. 
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The  behavior  of  the  imperfect  plate  (e  +  0)  is  best  visualized  by 
using  the  auxiliary  variable  y  =  (8+e).  The  equation  for  X  becomes  then: 

d5X  =  d4-y2  +  3d7^  +  td5^cr  "  ^d4E  +  2d7e^E^  "  ^d5Xcr  +  d7ee'>Y 

The  additional  effects  of  the  imperfection  are  seen  to  be  as  follows: 

(i)  A  hyperbola  which  is  mainly  proportional  to  e  is  added  to  the 
bifurcation  parabola. 

(ii)  The  A  coordinate  axis  shifts  £  units  to  the  right 

(iii) The  parabola  translates  downward  by  an  amount  which  is  essentially 
proportional  to  c 

Of  the  three,  the  last  one  is  the  least  important  and  may  be  omitted  for 
sufficiently  small  .  Indeed,  for  sufficiently  small  e  the  approximation 
which  preserves  all  the  important  features  is : 

X  -X^r^X  “  c/y)  +  -J-  (d^Y2  +  3d^ey) 

Figure  3  presents  the  equilibrium  paths  for  a  plate  with  the  following 
properties  (corresponding  to  a  plate  for  which  test  data  is  available) : 


y  =  0.5;  A  =  0.179;  I  =  0.694;  e  =  0.2 


For  this  plate:  A  =  4.281.  Also: 
cr 


A’  -  9.408 
cr 


A  ,  -  8.616 

min 


The  equilibrium  paths  are  shown  in  the  figure  for  the  following  im¬ 


perfections:  e  ■  0.2(4>q  =  -0.51);  e  *  0 (4>q  ■  -0.71);  e  ■  -0.01($o  ■  -0.72); 
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e  =  — 0 - 1  ( ct>o  =  -0.81).  The  bifurcation-type  behavior  for  e  =  0  and  the  post- 
buckling  unstable  '  chavior  for  e  =  -0.01  are  clearly  seen.  For  the  other 
imperfections  the  equilibrium  paths  are  stable  over  the  entire  range. 

Several  annular  plates  having  similar  configurations  were  tested 
in  the  past  few  years  [8,10,11].  The  results  of  the  tests  show  a  be¬ 
havior  which  is  very  similar  to  theresults  given  here.  In  addition, 
in  some  of  the  tests  jump  phenomenae  occurred  with  sudden  changes  in  the 
measured  rotations  and  subsequent  stable  behavior.  This  jump  is  very 
well  modeled  by  the  interaction  analysis  presented  here. 

10-  SUMMARY  OF  RESULTS 

(1)  The  capabilities  and  potentialities  of  semi-inextensional  plate 
and  shell  models  were  demonstrated  in  the  following: 

(a)  Development  of  the  ordinary  differential  equation  for 

the  semi-inextensional  axisymmetric  deformations  of  shells  of  revolution. 

(b)  Comparison  of  an  approximate  semi-inextensional  model  for 
the  annular  plate  with  a  corresponding  (but  more  complex)  extensional 
model. 

(2)  A  study  of  the  interactions  between  stiffener  eccentricities  and 
surface  imperfections  was  initiated  using  simple  models  of  the  annular 
plate. 

(a)  The  appearance  of  unstable  postbuckling  behavior  and  sen¬ 
sitivity  was  noted  and  examined.  The  sensitivity  and  its  corres¬ 
ponding  jump  were  restricted  to  a  selected  range  of  the  parameters 
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and  was  limited  in  magnitude.  It  is  possible  that  many  "small  scale" 
jump  phenomenae  which  occur  in  tested  stiffened  structures  are  of 
this  type.  As  demonstrated  here,  these  need  not  be  dangerous  but 
their  occurrence  should  be  taken  into  consideration. 

(b)  The  basic  formulas  for  a  more  detailed  study  using  a  more 
complex  model  were  developed  and  presented  in  a  manner  which  is 
suitable  for  analysis  (see  Appendix). 
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APPENDIX:  Equations  for  a  Three-Parameter  Semi-Inextens ional  Model 


By  introducing  an  additional  rotational  degree  of  freedom  into  the 
semi-inextensional  model,  its  capabilities  can  be  greatly  increased  so 
that  even  plat  es  with  smaller  central  openings  (0  •  ^  •  1)  or  heavier 

rotational  edge  restraints  (i  -  Oil)  and  higher)  can  be  successfully  treated 
and  greater  accuracy  can  lie  achieved.  fhe  penalty  is  in  the  increased 
amount  of  algebraic  complex  i.ty  The  basic  equations  are  developed  here 
and  brought  tc  a  form  which  con  be  treated  by  simple  algebraic  means. 

The.  question  as  to  whet  net  it  is  worthwhile  to  mvestigare  this 
approach  further,  or  Co  utilize  existing  nonlinear  programs  (such  as 
B0S0R  5)  instead,  is  Jett  open. 

The  model  is  def  ined  by  the  two  angular  parameters  in  the 

expression  for  the  rotation 

ifl  -  j> ,  +  p_( t /a) 

l  i 

and  by  the  displacement  u  cl  the  inner  edge  of  the  plate.  Using  (29), 

(7),  (8b),  the  expressions  tor  the  strain  measures  are: 


t.  “  0 

r 


u. 


K  <=  a 


Ke  '  +  7(W 


-34- 


Introducing  into  equation  (31)  for  the  total  potential  and  performing 
the  required  integrations,  the  following  expression  for  the  total  potential 
V  results: 

iv  -  a2C<2jJ<l-«6>*5  +ij<lV)»!<5+2-2^)  +  UV>  [£Vl(*l'+o)- 

-  Yg4>2w]  +  -|(1-P2)  [(4>2“<J>q)2  -  2$1<t>2w]  -  -|(l-y)  (4>^-<!>2)w  -  w22.np}  + 

+  D[  (1+v)  (1-p2) if)2  +  2(l+v)(l-p)<t>2(,j)i-(j)o)  -  (<|)1-$o)2Jlnp]  + 

+  -g-E{l(p<j>2+<j)1-<?o)  c  +  [ub  +  ( p<f>2-J-<|)1— d=0) e ] 2 A >  - 

-  a2P[^2  +<(>2-<j)2-w]  (A. 2) 

where: 

uh  1 

w  =  2 —  +  -rp3<t>2  +  +  u(<fi?-<t>2) 

a  J  2  12  1-0 

The  total  potential  can  be  now  minimized  with  respect  to  the  deformation 
parameters  (<J>2 ;  ^  ;w) .  The  results  of  this  process  lead,  after  some  mani¬ 
pulations  to  the  following  three  equations  for  the  parameters  in  terms 
of  the  loading  parameter  *  and  the  geometry: 

(|A5-4A2A3m)^  +  (|A4-18A2m-16A1A3m)^?1  +  (2A;3-36A1A2m)?2(?|-^) 

A2 

+  (12m3A-72A  A„m)4>07?  +  (6A„-144Ai;mH  (Tr-72)  +  (6 — Pirt+16A,m-l6y 3)e<j>^ 

1221  2  lllo  —  J  2 

A 

2A1Pm  _  A2m  _  _  _  _ 

+12  (— -  +  6A2m-3p2)e$2!t>1  +  6(~ - 2uz )  esf>2  (4>^— ^0>  +  12(A^m-u)e(^-^)  + 

A  A 

A-m  _  _  _  _ 

+24(~-  -p)e*1($1-?0)  +  2 [  (1+v)  (1-y)  +yIcJ*2  +  2(Ic-imi)  (^-^J  -  [(l-6A2m)$2 
A 

+2(l-l^A1m)?1  +  24me] A 


(A.  3) 
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■^•<A6— SA^m)  +  ■|(A5-10A2A3m)i|<t>1+  (A^-lbA^A^m)  <J>2<?^-!^)  +  |■(A^-12A2m)4>2'i>^  + 

+  2(A3-18A1A2m)ij!1(^-(t>^)  +  [^^(6+^)-  6y3]uei"|  +  6[A  m(12+  ^)— 4vi2 ) ye4>24>1  + 

A  A 

mA-  _  _  _  mA-  _  _  _ 

+  8( — —  -P3)e4>, (*.-*)  +  12(12mA.-v)ue(^-P)  +  6(—  -2p2)e<k  (♦.-*)  + 

3A  210  1  10  A  110 

+  2[(l+v)(l-u2)4-p2TcH2  +  2l(l+v)(l-p)+  Pl*cl  (T1-Te)  =  [-|(l-4inA3)?2  + 
(l-mA0)<j>^  +  24pme]  A  (A.  4) 

and : 

-ub  =  mX  +  ■^<V,-4A3m)'i'2  +  -^C1  ■' --6A0m)  +  •|<P-12A.]ra)  +  ^^4>2  + 

A 

+  —  (*.-*>  (A.  5) 

A  1  0 

Hire : 

A  =  1  -  p'(l-jA)  (i  1 - 6) 

ir.  =  Jj(A  -  £np)  ^ 

Also,  I  is  the  modified  moment  of  inertia  of  the  r  ing  with  r  espect  to 
the  midsurface  of  the  plate,  given  by: 

L  =  I  +  Ae?(i-m) 

c 

Other  notations  are  as  in  the  previous  chapters.  For  example: 

_  El 

la  =  bD1  *2  "  *2  a/h  etC’ 
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Kquations  (A.  3)  and  (A. 4)  have  been  derived  by  first  t«aking  V,^ 
and  V  ,,  respectively,  and  then  eliminating  by  using  (A.  5)  which  has 
been  derived  ( rom  V,w.  In  this  fashion  the  rotational  parameters  are 
decoupled  from  and  can  lie  solved  lor  separately.  The  general  form 
ol  Lhe  equal  tons  is: 

( i  l,l) 

wiu're  I.  and  M  arc  cubic  and  linear  operators,  respectively,  in  '{„  and  <K  . 

it  l  L 

IKiue,  ii  is  possible  to  eliminate  \  iron)  the  equations  and  arrive,  at  the 
s  ;  nglt'  eqra  i  i  on  : 

OlT.iTj)  -  U.,M1  -  i.^uT.,;^)  =  0  (A. 6) 

which  describes  the  internal  relationship  between  and  f  ,  irrespective 
of  the  li*ti ill. i,.  it  is  a  quaitic  equation  in 

I' .-•mg  assumed  values  lot  *  (starting  i  rom  =;^),  the  corresponding 
solutions  t  or  >,  c,;n  be  found  from  (A. 6)  which  becomes  a  fourth  order 
“olyiu'mtal  i n  the  f^(the  luadanientai  solution  is  the  one  starting  from 
-  0,  but  the  other  -elutions  can  be  useful  tor  exploring  additional 
equilibrium  h.anches).  for  any  pal  :  of  (1  the  loading  can  be  found 

from  (A. 1 ) .  T’iis  yields  the  equilibrium  path  which  can  be  represented  as 
=  (ifp)  where  (•  is  the  rotation  at  any  convenient  point  in  the  plate. 

The  only  mathematical  tool  needed  is  the  solution  of  a  fourth  older 


polvnomial  equation. 
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Figuce  2:  Deformation  and  Equilibrium. 


